AMENDMENTS TO THE SPECIFICATION 

Please amend the following paragraphs: [0011], [0046], [0057] and [0059]. No new 
matters have been introduced. 

[0011] The present invention discloses a system, method and software product to 
numerically simulate highly compressible material such as foam with a pragmatic 
approach. Under the assumption of uniaxial loading and isotropic material, a 
method of calculating stress function ^X) is developed as shown is FIGS. 2 and 3. 
The method eliminates the requirement of fitting a polynomial function; instead is 
calculated using only a set of the stress-strain curve from a uniaxial loading test. 
Based upon the data of this stress-strain curve, the stress function f(?i) is pre- 
computed for a range of equally spaced stretch ratios X. These results are stored in 
a computer's RAM as a lookup table and retrieved later during the solution phase of 
FEA. The present invention does not require time-consuming trial-and-error process 
of fitting polynomial coefficients. This is a huge advantage over the •^--^ \ ^ \ 
method. The stresses corresponding to a particular strain or stretch ratio can be 
interpolated rapidly from the pre-computed table directly during the lengthy solution 
phase. 

[0046] Elastic material is a material that. recovers to.its original shape after 
unloading (i.e., the load is removed). 

[0057] By rearranging formulas 240, 250, 260 and the rest of higher order 
formulas, the resulting equation for stress function 1{X) is rewritten as listed in 

3ilCI310. The equation 3-2#-310 js rewritten as a '\ o t form equation as shown in 
^3§320. Wherein the value of stress function f(X) for the stretch ratio of interest X 
equals to the summation of a sequence of a.^-"^ fT„|i'-'"' ) , where j is an integer related 

to the j-th term of the sequence, v is Poisson's ratio of the compressible material, 
and cr„(/i' is the stress value at stretch ratio i""' . Using the relationship between 
strain e and stretch ratio X in 340330. the equation 33Q-320 transforms to equation 
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S§0340. which can be evaluated solely with the stress values corresponding to the 
strains from the strain-stress curve obtained via a uniaxial tension/compression test 
of the material of interest. Finally, the nominal stress go and true stresses a are 
rewritten and dependent only on the stress function i{X) listed in formulas 350 and 
360, respectively 

[0059] With reference now to FIG. 5, it shows a flow chart 500 for computing the 
infinite series in equation Because the order of magnitude decreases 

drastically from one term to the next, the sequence converges very rapidly. 
According to one embodiment of the present invention, the flow chart 500 of 
evaluation of function f(A,) is summarized in FIG. 5. At 510, f(X) is assigned a value 
equals to X multiplied by ao(>--1) for a given X. The stress ao value at strain so or X-1 
is from engineering test strain-stress data (e.g., the curve defined in the user input 
phase of FEA software). At 520, A, is stored into a variable Xoid- A new variable Anew 
is stored as at 530, where v is Poisson's ratio of the compressible material. At 
540, a comparing test is performed for the absolute value of (Xnew-1) being less than 
or equal to a threshold value ^. In one embodiment, ^ is set to 0.01 . If the test 
succeeds, the computation of f(X) has finished, the rest of the terms in the infinite 
series is too small to affect the final result of the computation. If the test fails, 
another stress value at strain ?inew-1 is multiplied by Xnewand accumulated into 
function f(A,) at 550. At 560, the value of >.new is stored into Xoid- The process goes 
back to 530 until the computation finishes. 
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